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ABSTRACT 

Feedback from active galactic nuclei (AGN) has become a major component in sim- 
ulations of galaxy evolution, in particular for massive galaxies. AGN jets have been 
shown to provide a large amount of energy and are capable of quenching cooHng flows. 
Their impact on the host galaxy, however, is still not understood. Subgrid models of 
AGN activity in a galaxy evolution context so far have been mostly focused on the 
quenching of star formation. To shed more light on the actual physics of the "radio 
mode" part of AGN activity, we have performed simulations of the interaction of a 
powerful AGN jet with the massive gaseous disc (10^^ of a high-redshift galaxy. 
We spatially resolve both the jet and the clumpy, multi-phase interstellar medium 
(ISM) and include an explicit star formation model in the simulation. Following the 
system over more than 10^ years, we find that the jet activity excavates the central 
region, but overall causes a significant change to the shape of the density probability 
distribution function and hence the star formation rate due to the formation of a blast 
wave with strong compression and cooling in the ISM. This results in a ring- or disc- 
shaped population of young stars. At later times, the increase in star formation rate 
also occurs in the disc regions further out since the jet cocoon pressurizes the ISM. 
The total mass of the additionally formed stars may be up to 10^^ M© for one duty 
cycle. We discuss the details of this jet-induced star formation (positive feedback) and 
its potential consequences for galaxy evolution and observable signatures. 
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1 INTRODUCTION 

Large scale cosmological simulations often include AGN 
feedback in massive galaxies in order to avoid both fur- 
ther growth of massive galaxies via star formation and 
the presence of a young stellar population at low redshift 
(e.g. ICroton et all 1 20061: iBower et al.1 l20Q6l : iHopkins et"all 
I2OO6I : ISchawinski et al.ll2006l ). This approach is quite natu- 
ral since AGN activity is frequently observed in these galax- 
ies ([Best et al.|[2Q0 5V So far, the models focus on negative 
feedback (suppression of star formation) to quench cooling 
flows, remove potentially star- forming gas by expelling it 
from the galaxy, or heating the interstellar gas to high tem- 
peratures in the halo to prevent new stars from forming. 
In contrast to radiation-based "quasar mode" feedback at 
high accretion rates, "radio mode" feedback by jets is pow- 



* E-mail: v.gaibler@uni-heidelberg.de 



ered mechanically by the kinetic ene rgy of the jet beam 
(IScheuerlll974l : iBlandford k Reeslll974l ) . For the case of cool- 
ing flows in clusters of galaxies and giant ellipticals, there 
also is meanwhile considerable evidence from observations 
of X-ray cavities that the energies deposited by the jets 
suffice in many cases to suppress the formation of cool- 
ing flow s in the ho t circum-galactic and intra-clu ster me- 
dia (e .g. lBirzan et al . 2004; Dunn et al. 2005; Raff ertv et al 
2006 ) . Furthermore, theoretical studies (e.^. IZanni et al 



20051 : iGaibler et al.ll2009l : IO'Neill Jonesll2010l ) show a high 
thermalization efficiency for this interaction. 

In this paper, we examine the impact of jet feedback on 
the star formation in massive, gas-rich galaxies at high red- 
shift (2; 2 — 3). It is not obvious why jet feedback, which 
acts in a negative manner on the hot and dilute circum- 
galactic X-ray gas, should also act similarly on the inter- 
stellar medium (ISM), since this is a multi-phase medium 
with much higher densities and a pronounced clumpy den- 
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sity structure in dynamic equili brium between ^as at vari- 
ous d ensities and temperatures Avillez Breitschwerdtl 
An alternative option is the model of positive feed- 
back, where jet activi ty actuall y triggers or enhances star 
formation in a galaxy (|Silkll2005l ). This is, however, generally 
applied on the galaxy scales and without explicit hydrody- 
namical modelling in large-scale cosmological simulations. 

Observational studies on this issue so far are inconclu- 
sive. Locally, one third of the radio galaxies show bright 
emission from warm molecular ^sls and very low to moderate 
star formation rates (I Ogle et al.l l2007l . l2010l : [Nesvadba et al.l 
I2OIOI ). The jets seem to heat the cold molecular gas and 
keep it warm during the lifetime of the jet although the 
cooling times are much shorter than the jet lifetime - po- 
tentially due to the action of shocks, turbulent dissipation or 
cosmic rays (Ferland et III I2OO8I . I2OO9I : iGuillard et al. I I2OO9I : 
IPapa dopoulos et al. 2010). The jets in those sources may ac- 
tually induce formation of molecular gas b ut not cause the 
formation of a detectable number of stars (|Nesvadba et al.l 
l201lh . Shock- heated, dense and turbulent gas has also been 
predic ted in multi-ph ase turbulence simulations for jet co- 
coons (|Krause Alexander 2007i ). Positive feedback, on the 
other hand, is directly support ed by only a f ew ob servations 
locally ( Minkowski's o b ject i n Croft et al. (I2OO6I ). Centau- 
rus A in iMould et all feoOOl): iMorgantil (I2OIOI )). at inter- 
mediate fPKS225 0-41 inllnsk p et al l (l2008|~nd h i^h red- 
shift (4C 41.17 in lDev et al.l (|l997l l: iBickneh et al.l (|200Q \V 
But several studies point towards a link between jets and 
star formation, such as a high incidence o f recent star 
formation activity in compact radio sources (|Dicken et al.l 
I2OI2I ). the detection of large masses of cold molecular g,d .s 
in high redshift radio galaxies (e.g. lEmonts et al.l l201lh. 
alignment of excited CO with iets (iKlamer et al. 2004), 
alignment of continu um emission (|McCarthv et al.l il987l : 
IChambers et al.| 1987) and a correlation o f emission line 
gas ([McCarthv. van Breugel Kapahilll99lh with the radio 
source. Furthermore, several authors find evidence for ^ 



source, l^urthermore, several authors tind evidence tor young 
stellar population s in radio g:alaxies (iTadhunter et al" I2OO2I: 



Wills et al.1 I2OO2I : iBaldi Cai^^ '2008: 'Tadhunte r et al.l 
201lh and^heir central regions (Aretxaga et al. 2 00 y, pos- 



sibly accounting for a substantial mass fraction (|Holt et al.l 
I2OO7I ). It is important to keep in mind that the circumstances 
of jet feedback vary largely throughout these studies since 
gas masses and the velocity dispersions are much larger at 
high redshift (a setting that we are addressing in this paper) 
than in the local Universe, and variations amongst radio 
galaxies of a certain redshift regime are large, too. 

We are still in an early stage of understanding jet 
feedback, despite very detailed observational studies in 
this field. The lack of physical modelling and understand- 
ing of jet feedback is currently one of the major short- 
comings of "radio mode" feedback scenarios on galaxy 
scales. Only recently, high resolution hydrodynamical simu- 
lations in cluding a clumpy interstellar medium have become 
feasible (ISutherland Bicknelll l2007l : ICaibler et all I2OIII : 
I Wagner Bicknell[201l V Indeed, these studies have shown 
that the multi-phase structure of the ISM causes the mor- 
phology and jet-gas interaction to be qualitatively differ- 
ent from simulations without the cool er and denser phases . 
lAntonuccio-Delogu Sil^ (l2008h and iTortora et al.l (l2009l ) 
have studied the interaction of a powerful jet in 2D slab 
geometry with cool clumps of gas dispersed within a hot 



medium, finding only a slight increase in the star forma- 
tion (10 - 20 per cent) followed by a much stronger de- 
crease (more than 50 per cent) within a few million years 
due to heating and ablation of the cloud gas. However, the 
2D approach results in a different pressure evolution of the 
coc oon and a larger eff e ctive surface of the clouds than in 
3D. IWagner Bicknelll tOllh studied the interaction of a 
relativist ic jet with a clumpy gas distribution in a galaxy's 
centre at high resolution. They found a considerable impact 
of the filling factor of the dense gas on the interaction effi- 
ciency and that the dispersal of gas at large velocities may 
actually inhibit star formation. Yet, their simulated time- 
scales of < 1 Myr allow only a weak impact of cooling and 
the simulations do not explicitly include a star formation 
model. 

In our pr evious p aper, ([Caibler, Khochfar Krausd 

I2OIIL hereafter IPaper ih . we found that the action of jets 
in a clumpy medium results in a complex interaction with 
both individual clumps and the average density distribution 
of the gaseous distribution and results in stochastic asymme- 
tries of the jets as one signature of this strong interaction. 
The resulting blast wave in a gas-rich galaxy compresses 
structures to high density where cooling is highly effective 
and impedes the cloud's pressure support and thereby its 
adiabatic re-expansion. In this paper we extend our earlier 
study to include star formation. We will here focus on the 
impact of the jet on the star forming gas in the galaxy and 
its observable signatures. We find that a strong increase in 
star formation rate (SFR) is expected with a clear morpho- 
logical signature in the stellar distribution. While the total 
number of stars in a given galaxy probably is not changed 
much by this, it would still be expected to be observable 
while it is ongoing and can have a significant impact on the 
galaxy's evolution. 

The paper is organized as follows: in Sect. [2] we summa- 
rize the setu p of the simulation and state the differences to 
the setup of IPapei^ . The results are given in Sect. [3l sep- 
arating three different stages found in the simulation. The 
results are then discussed in Sect. 2] 



2 THE MODEL 

We model the host galaxy as clumpy, disc-like distribution of 
gas with a log-normal probability distribution of the density 
(median lOmpCm"^, mean 16.5mpcm~^), an exponential 
radial profile with scale radius Ro = 5 kpc and a sech^ ver- 
tical profile with scale height ho = 1.5 kpc, reminiscent of 
a massive and thic k galactic gaseous disc at high redshift 
(|Genzel et al.ll2006l V The width of the log- normal density 
probability distribution corresponds to a turbulence of Mach 
number 5 (^ 80 km s~ ^), which i s similar to large evolving 
discs at high redshift (|Forster S chreiber et al. 2009). The 
two-point structure of the density field is described by a 
Fourier spectrum that follows E{k) oc k~^^^ for large wave 
numbers /c, but is damped towards smaller wave numbers to 
suppress inhomogeneities on scales larger than the disc scale 
height. The total mass of the disc is 1.5 x 10^^ M©. The disc 
is embedded in a homogeneous hot atmosphere with a den- 
sity of 0.05 nip cm~^ and a temperature of 1.15 x 10^ K. The 
setup is almost identical to the one described in lPaper"^ : dif- 
ferences are noted in this section. The jet is implemented as 
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Figure 1. Density volume rendering of the central part of the domain (32 kpc box length) at t = 14 (left) and t = 22 Myr (right). Only 
the z > half is shown and densities close to the ambient X-ray gas are transparent to give a tomographic view. The colour bars show 
log p in units of nip cm~^. 




Figure 2. Density volume rendering of the high density gas in the disc (face-on view) at t = 14 (left) and t = 22 Myr (right). The box 
length is 32 kpc, lower density regions are transparent and the colour scale is the same as in Fig. [T] 



cylindrical orifice of jet plasma (density 5 x 10~^mpcm~^, 
speed 0.8 c) with the same pressure as the environment, a 
radius of rj =0.4 kpc and an initial length in both directions 
of 3rj , respectively. 

We include radiative cooling by a tabulated cool- 
in^ function for atomic p rocesses down to 10^ K as 
[Sutherland Dopital (|l993l i with a metallicity Z = 0.5 Zq 
(lErbl 120081 )." which is a very important effect in the disc 
due to the high densities and correspondingly short cooling 
times. No strong dependence on the metallicity or the exact 
form of the cooling function is expected since the mechanism 
found in our simulations can simply be understood by the 
fact that the cooling becomes more efficient (shorter cooling 
times) for denser gas. To avoid the need to include a fine- 
tuned feedback recipe to stabilize the cooling disc against the 
vertical gravity component we exclude gravity in the simu- 
lations and stabilize the disc by imposing a minimum tem- 
perature of 10^ K, which mimics a strong radiative heating 
from the stars that prevents any cooling below the thresh- 



old although this might still be expected in dense, shielded 
regions. 

Since cooling is included, no truly static setup of the 
disc is possible and we run control simulations of the disc 
without the jet for comparison in order to mon itor the evo- 
lution of the undisturbed disc. In contrast to IPaper j . we 
use constant pressure for the initial conditions of the en- 
tire domain, however enforcing a pressure corresponding to 
a minimum temperature of Tmin//i — 10^ K where neces- 
sary (/i: mean particle mass in proton masses). Although 
the disc gas cools rapidly to Tmin and drops out of pressure 
balance, it shows only little evolution in the fluid variables 
with little impact on the dynami cs (howe ver, resulting in 
a different jet asymmetry than in IPape^^ . To allow for a 
more relaxed disc state, the jet of power Lkin = 5.5 x 10^^ 
erg s~^ is only started at t = 10 Myr. T he simulations are 
carried out using the RAMSES 3.0 code (lTevssierll20"oi ). a 
non-relativistic second-order Godunov-type shock-capturing 
adaptive mesh refinement code. The grid is refined to the 
maximum of 62.5 pc cell size in all regions of interest. For the 
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many dense clumps, the cooling length cannot be resolved 
and hence shells formed by shocks are expected to be thin- 
ner and denser in nature for these cases. We include tracer 
fields for both the jet plasma and the cool disc gas to fol- 
low its spatial distribution independent of mixing processes. 
Magnetic fields are not considered in this study. Our ear - 
lier magnetohydrodynamic simulations (|Gaibler et al.|[2009l ) 
showed that for the morphology of the jet cocoon at large 
scales they play an important role since they efficiently sup- 
press Kelvin-Helmholtz instabilities at the contact discon- 
tinuity. Although the interaction of the jet with individual 
clumps of gas may be affected by magnetic fields and reduce 
the efficiency of ablation, this is beyond the scope of this 
paper and we leave it to a future study. 

In this study, we use a star formation model 
(jRasera L Tevssie j [2OO6I ) . where particles are generated in 
higher density regions of hydrogen number density nu > 
{p > 1.31 mpTi^) according to = ep/t^^ where tff is 
the local gas free-fall time and e the star formation effi- 
ciency. This r nethod is widely used in hydrodynamic sim- 
ulations (e.^. [Kravtsovl l2003l: ISpr ingel Hernauist '2003|; 
ISchave Dalla Vecchial I2OO8I: IPubois Tevssier 200j), 
supported by observations ([Krumholz Tanll2007h and is 
able to recover the Kennicutt-Schmidt law ([Kennicuttll 19981 ) 
in simulations. After spawning a star particle, the corre- 
sponding gas is removed from the grid. We use two pairs 
of parameters (n^ = 0.1cm~^,e = 0.025) and (n^ = 
5cm~"^, e = 0.05) to check the impact of star formation be- 
ing more strongly localized in high density regions. The effi- 
ciency parameters e were chosen to get a SFR ^ 150 M© 
yr~^ sim ilar to those observed in massive high-redshift 
galaxies (jGenzel et al.l I2OIOI ). The higher star formation 
threshold can be regarded as being more realistic for our 
resolution, but shows more relaxation in the disc (which it- 
self tends to decrease the SFR with time). 



3 RESULTS 

We here briefly summarize the formation of the blast wave 
as reported in I Paper j . The jets produce a blast wave in 
the central region of the disc since the density contrast be- 
tween the jet and the dense gas (77 = Pj/pdisc) is very strong 
(?7 < 10~^), causing the jet head to stall and deposit its 
energy in mostly thermal form into the disc gas. This is 
also true at later times, when the jet breaks out of the 
disc, but since the density contrast is weaker, a stronger 
forward propagation results. Figs. [T] and [2] give an overview 
of the evo lution of the jet-disc interaction (see also fig. 3 
of IPape^ and the supplementary movie there). Since the 
jet fe edback is energ!:y- driven, the large mechanical advan- 
tage (jWagner Bicknelll201lh allows a much stronger effect 
then expected for a mo mentum- driven jet feedback model 
(iKrause k, GaiblerllioTol ). The early evolution hence is dom- 
inated by the effects of the blast wave rather than direct 
interaction of the jet beam with the ISM clouds, but the 
dense clumps inhibit the propagation of the jet and thus 
help to deposit the jet energy closer to the centre of the 
galaxy, thereby feeding the blast wave. 

The compression of gas by the blast wave and its short 
cooling times result in a considerable increase in the SFR 
despite the excavation in the centre (Fig. [3]). The simula- 
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Figure 3. Total star formation rate for star formation thresholds 
of = 0.1 cm~^ {top) and = 5cm~^ (bottom). The jet is 
started at f = 10 Myr. 



tion with a higher star formation threshold (n^ = 5 cm~^) 
shows a much stronger increase in star formation than the 
one with a lower threshold (48OM0yr~^ vs. 28OM0yr~^ at 
t = 22 Myr). This can be understood as a consequence of 
the stronger cooling and compression at these potential sites 
of star formation. Also, the star formation is more concen- 
trated towards the centre with the higher threshold due to 
the radial density profile. We have performed control simu- 
lations without the jets to ensure that the changes in star 
formation are actually caused by the jet activity, since the 
undisturbed disc cannot be in hydrostatic equilibrium due 
to the cooling. Indeed, the undisturbed disc evolves much 
more moderately than the disc with active jets. The relax- 
ation effects in the disc are somewhat larger for the higher 
threshold run since stars are formed in smaller regions with 
higher densities, but they are generally opposite to the ef- 
fects caused by the jet (they tend to lower the SFR). By 
the end of the simulation, the SFR is still increasing contin- 
uously. This can be understood by looking at the star for- 
mation in three different cylindrical regions around the disc 
axis (Fig. 2]). In addition. Fig. [5] shows face-on and edge-on 
maps of the stars at three different times. 

The SFR in the inner region (R < 3 kpc) drops after an 
initial increase caused by compression and the subsequent 
removal of gas. In the intermediate region, between 3 and 6 
kpc, the SFR increases strongly by a factor of 2.5 and 3.5 
for the low and high threshold, respectively, and then settles 
to this value since the main interaction region for the jet is 
then outside the disc and the increased SFR only continues 
eating up the generated high-density gas. The outer region 



© ??? RAS, MNRAS QOO.fUfM] 



Jet-induced star formation in gas-rich galaxies 5 



150 



100 



50 



[0, 3] kpc 

— — [3, 6] kpc 

[6, 16] kpc 



= 0.1 



^ — — 




10 15 
t [Myr] 



[0, 3] kpc 
[3, 6] kpc 
[6, 16] kpc 




[Myr] 



Figure 4. Evolution of the star formation rate for three cyhndri- 
cal regions of the disc: the central area (R < 3 kpc), a ring-like re- 
gion (3 kpc ^ R < 6 kpc) and the outer disc (6 kpc ^ R < 16 kpc). 
The top panel is for a star formation threshold = 0.1 cm~"^, 
the bottom panel for = Scm""^. Only stars with Ax < 10 
kpc from the disc plane are considered, the disc-only simulation 
is shown by blue, thin lines, the red, thick lines include the jet. 



(R ^ 6 kpc) only shows an increasing SFR at t > 14 Myr, 
but it continues to rise until the end of the simulation. The 
increase in SFR at large radii is generated by the overpres- 
sure in the jet's bow shock and cocoon, which is compressing 
the entire ISM of the disc (Fig. [6)) and thereby pushes its 
density beyond the star formation threshold. The control 
simulations show a slightly decreasing SFR except for the 
outer area, where it increases a bit since the ambient gas 
is already at slightly higher pressures than the cooling disc 
gas. So far, we cannot tell from the simulations how the long- 
term evolution of the SFR in the outer region is changed by 
the jet and how this depends on the modelling of the disc. 
Yet it is evident from the simulations, that this pressurizing 
of the disc occurs and may result in a large number of stars 
being formed. 

For the overall impact of the jet activity, we identify 
three phases that explain the evolution of star formation in 
the disc: (a) excavation of a central cavity, (b) compression 
of disc gas and triggering of ring-like star formation by the 
blast wave and (c) late-time evolution of the disc in the jet 
cocoon. In the following subsections, we will address these 
phases in more detail. The phases should not be understood 
as consecutive phases, since they overlap temporally. 



3.1 Phase A: Formation of a central cavity 

The most direct consequence of the blast wave is the forma- 
tion of a cavity in the disc centre since the bow shock pushes 
intermediate density gas outwards. The expansion can be ap- 
proximated by a spherical shape, and assuming a constant 
energy input by the jet, the blas t wave evoluti on can then 
be described by the one given in iKrausd (l2005l ). keeping in 
mind the expectable deviations due to the inhomogeneities 
in the disc. The blast wave expands until it leaves the disc 
vertically, then venting its high pressure efficiently though 
this channel since the weaker density contrast there results 
in a faster jet propagation. Consequently, the lateral expan- 
sion of the cocoon into the disc slows down. As can be seen 
from the simulation, this breakout happens approximately 
when the bow shock reaches a radius of rcav 2/io, where 
the vertical density profile dropped to 0.07 of the midplane 
density. The time needed for this is, according to the blast 
wave approximation. 



(2/io 



.5/3 



5Lk 



47Tp 



-1/3 



(1) 



with tcav ~ 3 Myr for our simulation parameters and ap- 
plying the mean central disc density as p. This means that 
for thinner discs and lower densities, the time-scales can be 
much shorter than 1 Myr, with rcav <C 1 kpc. 

The "cavity" actually is not empty but still contains 
a considerable amount of gas (Fig. (Tj): about half the mass 
initially present in the centre is still there after the passage 
of the blast wave. However, the gas is located in very dense 
clumps and filaments with only a small filling factor. 

We can furthermore estimate the energy deposited in 
this blast wave phase into removal of gas by assuming that 
half the mass of this volume is accelerated to a speed of 
V = 0.5rcav/tcav. With M ^ 167ihop/3 , we get 



^ki] 



(2) 



which is in our case gives 5 x 10^^ erg. For comparison: 



the jet injects 5 x 10 erg within tc 



: 3 Myr. 



In our simulation, we start with an initial jets length of 
already 1.2 kpc, so we actually have two blast waves orig- 
inating from the tip of the jet beam, respectively. These 
blast waves meet near the disc midplane, get shocked by 
their respective counterparts and compre ss ^as in to a thin 
disc-like structure (easily visible in fig. 3 of lPaper"^ , which is 
immediately destroyed by the shear flows in the cocoon and 
fragments into the filaments and clumps that remain located 
in the cavity. If the imposed initial jet length were zero, the 
two blast waves would be merged from the beginning as a 
single blast wave originating from the disc centre, skipping 
the intermediate creation of a dense, double-shocked cavity 
disc. We can now estimate what can be expected in nature 
for this early blast wave formation. 

The blast wave behaviour is a result of the jets hav- 
ing a much lower density than the ambient gas (77 ^ 1), 
where then momentum balance enforces a slow propagation 
and most of the energy is deposited in thermal energy. In 
contrast, high density jets (77 ^ 1) show a ballistic prop- 
agation behaviour. Hence, the blast wave phase will begin 
once the jet density drops strongly below the local ambient 
gas (i.e. once it encounters a region of much higher density). 
Larger jet beam widths are observed at larger distances from 
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Figure 5. Star formation maps for t = 10 Myr (top), t = 14 Myr (middle) and t = 22 Myr (bottom row) as RGB composite (left: 
edge-on view, right: face-on view). Newly formed stars (formed within the last 1 Myr) are shown in blue, all other are shown either in 
red (stars formed at t ^ 10 Myr, in the undisturbed disc) or green (stars formed at t > 10 Myr, jet was active). The images refer to the 
simulation with the high density threshold (n^ = 5cm~^), a scale bar with 5 kpc length is given. 
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Figure 6. Pressure slices through z = at t = 13 (left) and t = 22 Myr (right), showing logp normahzed to the ambient pressure. 
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Figure 7. Evolution of the mass inside a spherical region r < ro 
around the centre with ro = 2 kpc and 3 kpc. Initially overwriting 
the jet nozzle region accounts for 5.4 x IO^Mq, and < 10^ Mq over 
the entire run time. 



the galactic nucleus, stable jets showing only small opening 
angles of a few degrees. As jets expand slowly, their density 
decreases oc with distance (however, note that a Lorentz 
factor ^ 1 on the pc scale increases the effective density of 
the jet). Correspondingly, the jet might become underdense 
already at the scale of 10 pc to reach the strong density con- 
trast we have at the kpc scale. These numbers are clearly 
very dependent on the jet properties and the gas distribu- 
tion in the centre of the galaxy, but help putting this early 
stage into context. Once the two blast waves have merged, 
the evolution becomes independent of the initial blast wave 
locations. This consideration tells us that our setup is ro- 
bust and despite the expected morphological differences at 
small scales, the results of this study would not be affected 
significantly. 

Some disc gas gets expelled by the interaction with the 
jet, as can be determined by the compressible disc gas tracer 
field. A mass of the order of a few 10^ on each side (for 
|x| > 10 kpc) is moved outwards from the initial disc (|x| < 4 
kpc) with an asymmetric distribution (2 vs. 8 x 10^ M©), 
as expected for the asymmetric interaction with the disc 
gas. By t = 24 Myr, it extends out to ~ 20 kpc above 
the disc plane. We note that the star formation at larger 



distances from the disc (more than 10 kpc above the disc 
plane) is ignored in our SFR - it is generally low (lO'^M© 
for = 5cm~^ during the entire simulation). 



3.2 Phase B: Ring-like star formation 

While the blast wave creates a cavity in the disc centre, it 
pushes gas outward and compresses the disc at the cavity 
boundary with a ring-like or hour-glass shape. While the 
bow shock in vertical direction reaches out to lower den- 
sities, where the propagation becomes easier, it is much 
harder near the disc plane. The dense clumps in the disc 
with densities of > lOmpCm"^ have cooling times of < 100 
yr at temperatures slightly above T ^ Tmin- Once they get 
shocked by the blast wave, they are compressed but lose 
their pressure support almost instantly, whic h prevents later 
re-expansion (see also the simulations by iMellema et al.l 
I2OO2I : [Fragile et al1l2004l ). This happens even after the bow 
shock has passed since the post-shock pressure is ^ 10^ 
times higher than initially in the disc, causing the clumps 
to collapse further. 

The bow shock does not propagate very far into the 
disc due to the large mass in the disc and since the lat- 
eral expansion stalls once the jet has broken through the 
disc vertically. The strong compression of gas at the cav- 
ity boundary, however, leads to very efficient star formation 
with a ring-like morphology. Although our physical under- 
standing of star formation is still limited, it is clear that 
high gas densities and a strong increase in environmental 
pressure have a strong impact on the actual star fo rmation 
processes in molecular clouds (|Krumholz et al.lboilh , which 
can not be modelled within this simulation. 

The strong increase in density can also be seen clearly 
in the filling factor of the dense gas (Fig. [8|). Before the jet 
is launched, the filling factor for high-density gas decreases 
since this phase at T = Tmin has a higher pressure than 
the lower density gas and expands (this is the relaxation 
process). However, once the jet becomes active the filling 
factor increases by approximately one order of magnitude. 
The largely increased mass in dense regions then results in 
much larger star formation. 
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Figure 8. Evolution of the mass and volume filling factors of the 
disc. The time evolution is indicate by colour, going from blue 
via green and yellow to red. For early times (blue, dashed lines, 
t = 0, 5 and 10 Myr), when the jet is still off, the filling factor 
at p > lOmpCm""^) decreases. With an active jet (solid lines, 
t = 12, 14, 16, 18, 20 and 22 Myr) it quickly increases, but then 
decreases again for the highest densities. Please note the jump 
between t = 10 Myr and t = 12 Myr. The star formation threshold 
Ui, = 0.1 cm~^ is indicated as grey vertical line (dashed). The disc 
is defined as the geometrical region, where the disc density profile 
is larger than exp(— 2) = 13.5 per cent of its central value. 



The phase diagrams at three different times shown in 
Fig. [9] give a different view of the interaction of the jet with 
the dense and dilute disc phases. Two lines, indicating the 
path for adiabatic expansion and for constant pressure, are 
overplotted. Gas cooling to T = Tmin inevitably causes pres- 
sure gradients within the disc, resulting in the hot and the 
very dense phase of the ISM being at higher pressure than 
the intermediate density gas in the disc and consequently 
expanding during the early relaxation phase. Without ad- 
ditional heating processes (note that supernova feedback is 
not enough to balance the large radiative losses) , no dynam- 
ical equilibrium can be established to keep the disc stable. 
We deliberately have chosen to neglect these processes since 
we cannot model them in a physically correct way and our 
setup without gravity exhibits only weak relaxation during 
the simulated time. At t = 10 Myr, only little gas is left 
at intermediate temperatures. Once the jet is active (t = 12 
Myr), the jet cocoon appears as a new region at higher pres- 
sures. It extends from high temperatures, overpressured with 
respect to the jet nozzle and the ambient gas, down to lower 
temperatures and high densities. The latter part is the inter- 
action and mixing zone between the dense disc gas and the 
very low density jet plasma. Here, the cooling times are short 



and the gas cools down quickly. There is a clear accumula- 
tion of gas at T//J ^ 10^ K (where cooling times become 
somewhat larger) before it ends up quickly with T — ^min 
and underpressured, then reaching up to higher densities. At 
t = 20 Myr, the cocoon is less overpressured and visible as a 
distinct p — const feature. Now the interaction with the am- 
bient gas is more prominent and the gas distribution forms 
a downward-bent cooling path, which supplies the disc with 
an increased amount of dense, star- forming gas. 



3.3 Phase C: Later evolution in the jet cocoon 

At later times, when the jet mostly interacts with the am- 
bient gas, the changed environment of the gas disc becomes 
important: the bow shock now encloses the entire disc and 
the increased pressure of the jet cocoon affects the disc. This 
is due both to the ram pressure of the backflow with its vor- 
tices as well as the thermal pressure of the cocoon. Since the 
cocoon turbulence is in the subsonic or transonic regime, the 
thermal pressure dominates somewhat over the ram pres- 
sure. The cocoon pressure now sets the environmental pres- 
sure for the gas disc, which gets compressed. In our simu- 
lation, this is mostly on the surface of the disc, but with 
time this pressure will affect the entire disc. There is also 
ablation at the disc surface by the cocoon plasma, however 
this is mostly intermediate density gas - the high density, 
star- forming gas is much more robust and harder to affect, 
at sub-resolution scales it would also be stabilized by the 
self- gravity of the cloud. 

Because the cocoon pressure drops as the radio source 
expands , the disc compression effect wil l event ually come to 
an end. iGaibler, Krause k, Camenzindl (|2009h find that jet 
cocoons are not strongly overpressured anymore once the 
source extends over more than ^ 100 kpc, where self-similar 
expansion then breaks down. We expect disc compression to 
be effective over approximately the active phase of the source 
or a bit longer. The total stellar mass formed by phase C is 
comparable to that of the previously described phase B. 



3.4 Stellar velocities 

The dense gas is not only compressed by the blast wave 
and forms stars - it is also accelerated by the ram pressure, 
mostly in the radial direction. This motion is imprinted on 
the newly formed stars within the clumps, which thus can 
show rather high radial velocities fFig. llOp of several 100 km 
s~^, even exceeding the escape velocity of the galaxy. The 
disc-only simulations, in contrast, show much lower veloci- 
ties caused by the gas motions. This higher velocity disper- 
sion may actually be an observable signature of these stars. 
While we expect this to occur also in nature, we caution that 
the present simulations are not sufficient for the determina- 
tion of credible velocities because of two important points: 
first, the stars are not decelerated by gravitational forces 
since gravity is neglected in the setup, which is acceptable 
for the hydrodynamics on the short time-scales considered 
but has considerable impact on the stellar velocities. Second 
and even more important, stars form in the cores of molec- 
ular clouds at densities rather higher than the ones resolved 
here. Hence the real velocities of these cores are consider- 
ably lower than for the gas at densities corresponding to 



© ??? RAS, MNRAS QOQ.fTlfM] 



Jet-induced star formation in gas-rich galaxies 9 



-4-2 2 




og p [rrip/ cm'^] 

Figure 9. Phase diagram for the cells in the central (32kpc)^ box, which enclosed the entire disc: temperature Tf/j, vs. density p for 
three different times (the jet starts at t = 10 Myr. For comparison, lines of constant pressure as well as adiabatic expansion/compression 
(T oc p^~-^) are overplotted (dotted), the star formation threshold = 0.1 cm~^ is shown as vertical line (dashed) and the location of 
the jet nozzle is indicated. Mass-weighted density and temperature histograms with logarithmic bins are attached to the phase diagram, 
respectively. 
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Figure 10. Velocity distribution of the stars at t = 22 Myr for 

the low (rii, = 0.1 cm~^, top) and high (n^ = 5cm~^, bottom) 
star formation threshold. The contributions are shown separately 
for stars formed before (dotted) and after (solid) the jet starts at 
t = 10 Myr. The histograms show the mass of the stars within a 
bin size of 20 km s"-*^. 

the star formation threshold . This is expected to depend 
on the internal structure of the clouds and the efficiency 
of momentum transfer. Accordingly, the measured velocities 
are upper limits. We note, however, that the velocity distri- 
butions for the two star formation thresholds considered are 
not drastically different. About 8.6 per cent (n^ = 0.1 cm~^) 
or 8.3 per cent (n^ = 5 cm~^) of the stars (2 - 3 x 10^ M© in 
absolute numbers) formed at t ^ 10 Myr have velocities of 
^ 300 km s~^. Follow-up simulations with gravity and much 
larger resolution for the gas clouds will be needed to assess 
this aspect quantitatively. 



4 DISCUSSION 

4.1 Consequences for galaxy evolution 

With the present simulations, we examine the impact of jets 
on a massive disc-like galaxy. The properties of the gaseous 
disc were chosen to be comparable with (non-merger) mas- 
sive star- forming galaxies observed at redshifts of z > 1. 
The onset of the jets lead to an increase in SFR by a factor 
of two and more. During 12 Myr of jet activity, an addi- 
tional stellar mass of 1.6 X 10^ Mo (7.2 X 10^ M©) is formed 
for our runs with a star formation threshold of = 5 cm""^ 
(0.1 cm~^), respectively, compared to the quiescent disc sim- 
ulation. Since the simulations are covering only ^ 12 Myr 



of jet activity, it is still difficult to estimate the long-term 
effects on the galaxy. As for the generally assumed negative 
feedback from jets, we have found that the reduction of star 
formation occurs mostly in the central region of the disc 
and is much smaller than the enhancement of star forma- 
tion in the rest of the galaxy. By the end of the simulation, 
the SFR at intermediate radii does not increase further and 
would decrease slowly as the densest clumps are converted 
into stars. At larger radii, the SFR is still continuously in- 
creasing and is expected to stay at a high level until the jet 
cocoon pressure h as dropped to a valu e near the original am- 
bient pressure (cf . iGaibler et al]l2009l ) , which happens after 
a few 10^ yr. Altogether, the total stellar mass generated 
additionally by the jet activity may be up to ^ 10^° M©. 
While this is only a fraction of the total stars formed con- 
sidering the high quiescent (disc-only) SFR, these stars are 
distributed non-uniformly throughout the disc and are of 
similar age. 

Both SFR-enhancing phases (B and C) are caused by 
the large amount of thermal energy generated by the jet as it 
forces its way though the ISM. Although direct interaction 
of the jet beam with clouds is limited to a small volume, the 
resulting pressure affects a fair fraction of the disc at the 
early stage and eventually all of the galaxy once the bow 
shock has reached beyond the galaxy's radial extent. By 
that time (in our simulation within ^ 10 Myr), the galaxy 
is embedded in a hot and low-density atmosphere (the jet 
cocoon) with a pressure about one order or magnitude larger 
than initially. It is quite evident that increasing pressure in 
and around a galaxy within a relatively short time will have 
a considerable effect on its star formation, causing higher 
densities and higher star formation rates eventually - inde- 
pendent of whether star formation rates are understood in 
terms of the local free-fall time (as done by the code), the 
Kennicutt-Schmid t law ( Kennicutt 1998) or a global depen- 
denc e on pressure (IWong Blitjl2002l : lBlitz Rosolowskvl 
|200^ 

IWagner Bicknelll (l201lh note that the ablation time- 
scale in general would be shorter than the cloud collapse 
time-scale due to the large velocities in the jet cocoon, and 
destruction of the clouds hence would be expected without 
any possibly associated star formation. However, while abla- 
tion may actually remove a considerable fraction of the dense 
clouds, we expect that only the cloud envelopes with their 
lower densities would be affected. The star formation takes 
place only in the highest density regions where, at densities 
of ^ 10^ cm~^ and considering stabilization by self-gravity, 
the jet will hardly have a totally destructive effect, in the 
worst case leaving a naked high-density core. Shocks driven 
into the clouds will be strongly radiative and dense shells 
formed by that would additionally p rotect the core even 
thou gh the c loud may fragment (cf. iMellema et al. I I2OO2I : 
Coo per et al.|[2009i V Since the environmental pressure of the 
cloud jumps to very large values, the cloud core will eventu- 
ally be compressed to a higher density with an anticipated 
increase in star formation there. This might even happen in 
the cleared-out central disc region, where naked molecular 
cloud cores could be left embedded in the jet cocoon. This 
can be further quantified by the time scales of the processes 
competing with the cloud collapse. Following ICooper et al.l 
(2009) and using the data from the simulation and a fidu- 
cial radiative cloud of lOOmp cm~^ and 100 pc radius, the 
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cloud crushing time is > 10^ years and the growth time scale 
for Kelvin-Helmholtz instabilities is 10^ — 10^ years. These 
are considerably larger than the cooling time for the dense 
clouds, and the Kelvin-Helmholtz time scale would even be 
larger if magnetic fields were present in the clouds. Hence 
the cloud would be stable against ablation until self-gravity 
takes over with a free-fall time scale of 3 Myr and subse- 
quent star formation would be expected. 

If star-forming molecular clouds are seen as transient 
objects generated within a dynamic ISM of different phases, 
one could argue that destruction of warm (intermediate den- 
sity) gas by the jet as in the cloud envelopes cuts off the 
supply necessary for the formation of molecular clouds and 
thereby would lower the star formation rate on longer time- 
scales. In the present simulations, however, the amount of 
strongly cooling gas is not changed much by the jet activity: 
depending on the selection threshold, it decreases only by 
3 per cent for gas with cooling times Tc < 10^ yr or even 
increases by 10 per cent for Tc < 10^ yr between t = 10 and 
t = 20 Myr, hence providing a reservoir for molecular cloud 
formation also after the jet activity has ceased. 

We have performed the simulation with two different 



0.1 cm and 5 cm 



to 



star formation thresholds (r 
examine resolution effects that are clearly expected for a 
setup with strong cooling. For both thresholds, the local 
Jeans length is resolved by ^ 20 cells, and it could even be 
raised to 100 cm""^ for 4 cells per Jeans length. A stronger 
increase in SFR is found for the higher threshold, where 
star formation is more concentrated to the densest regions 
of the clouds, which we expect to be more realistic. Fur- 
thermore, the change in the density probability distribution 
function (PDF) caused by the onset of the jet activity is a 
robust and purely hydrodynam ical result (cf. compressional 
forcing, ^Federr ath et al.ll201Ql ) with the direct consequence 
of a larger expected star formation rate independent of the 
actual implementation; the time evolution of the PDF is 
much larger than the differences between the runs with dif- 
ferent thresholds at equal times. This increase in the den- 
sit y PDF at high dens i ties is also found in the simulations 
bv lWagner Bicknelll (|201lL private communication), with 
some differences due to the different setup. We are hence 
confident that the results are robust, and the enhancement 
of star formation effects might even be stronger if self-gravity 
of the gas clouds were considered. 

Thermal or kinetic AGN feedback in cosmological and 
galaxy evolution simulations (modelling radio and/or quasar 
mode) is often thought to heat or expel gas from the 
gala xy and to be acting mostly on the cir cum- galactic gas 
fe.g. ISpringel et al.|[2005l: ICroton et al. ' 2006; S iiacki et al 



20Q7l : lBower et al.ll2008l : lBooth Schava2009 - Dubois et al 



201l|), lowering the star formation rates. Our setup, in con- 
trast, has shown the entirely opposite effect: the formation 
of an additional population of stars. The need for an addi- 
tio nal boost of star formation has been recently suggested 
bv lKhochfarfcSilki d 201lh . They argue that to recover the 
high specific star formation rates that are observed in galax- 
ies at z > 5, ACDM models need to introduce stochastic 
boosts in star formation. Such boosts could indeed be trig- 
gered by jet activity as presented here in gas-rich galax- 
ies at 2; > 5. The main reason for the difference between 
our results and common AGN feedback simulations lies in 
the different treatment of the ISM and the collimated jets: 



whether it forms a rather continuous gas distribution or in- 
dividual small-scale clouds are resolved and interact with 
the propagating jets. While on the AGN activity time-scale, 
our simulation shows a considerable positive feedback, this 
does not yet imply that there could be no negative feedback 
associated with jets: First, it is important to keep in mind 
that outside the galaxy, in the rather uniformly distributed 
and diffuse circum-galactic or intra-cluster gas, the feedback 
is still negative, heating the gas and quenching cooling flows. 
Second, on the galaxy scale, the formation of ^ 10^° M© of 
stars (almost instantaneous compared to galaxy evolution 
time-scales and concentrated in certain regions, cf. Fig. [5]) 
will eventually be joined by the stellar feedback from these 
stars. This may have a strongly negative effect on the later 
evolution of star formation. For smaller gas masses, the pos- 
itive effect might be very small and more efficient ablation 
might cause negative feedback to be dominant for the lower- 
density environment, potentially resulting in positive feed- 
back in gas-rich (high-z) and negative feedback in gas-poor 
{\ow-z) environments. We note that while this study ad- 
dresses the physics of ongoing jet feedback, it is still unclear 
under which circumstances AGN jets are launched and how 
frequently jet activity would occur in gas-rich galaxies. 

Clearly the present simulations do not include all 
physics that might be relevant; magnetic fields and cosmic 
rays might play a role, and small-scale shocks or turbulent 
mixing layers may be effective below the finite resolution 
of our hydrodynamical simulation. However, these effects 
would affect discs without jets in a similar way unless the 
star formation mechanism is particularly sensitive to the 
changes. Assuming that giant molecular clouds are essen- 
tially self-regulated inside fe.g. iKrumholz et al.ll2009l V it is 
at least not obvious why the molecular cloud processes in 
the core would change drastically, except that the hydro- 
dynamical interaction with their environment changes the 
global cloud properties. In our simulations, compression by 
the blastwave and the jet cocoon in general is the dominant 
effect, with ablation and stirring of the cold clouds being 
limited to the surface of the disc. The long-term evolution 
of the galaxy and an exploration of the parameter space are 
beyond the scope of this paper and should be addressed in 
future studies. 



4.2 Comparison to observations 

As noted in Sect. [U the observational evidence for positive 
feedback or jet-induced star formation is still limited, while 
the scenario in the light of the presented simulations appears 
almost inevitable from a theoretical point of view. Clearly, 
gas masses and jet powers in the local Universe vary over a 
large range and the present study only probes one point in a 
large parameter space (powerful jets in massive high-redshift 
gaseous discs). Accordingly, the star formation induced by 
the jets in the often quite different nearby objects will be on 
a considerably smaller level and possibly be unrecognized. 
We have seen in the simulation that the jet causes a mor- 
phological disturbance in the stellar distribution - a ring or 
disc of stars is generated by the early blast wave of the jet 
expanding in the galaxy. Since the expansion becomes much 
slower and stalls once the jet is able to break out vertically 
through the dense gas, the radius of this stellar population 
should be on the order of a few vertical scale heights of 
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the gas. For galaxies in the local Universe, this would be 
at som e 100 pc. Interestingly, ! Jackson, Tadhunter Sparksl 
(| 19981 ) found a ring of young stars around the centre of 
Cygnus A, the most powerful radio source at low redshift 
(z — 0.056). They argue that the blue condensations are 
almost certainly due to young stars and not caused by scat- 
tering of the obscured quasar; the continuum energy distri- 
bution can be modelled by starburst activity much less than 
1 Gyr ago, possibly an instantaneous burst with t ^ 10 Myr, 
which is comp a tible with the radio source age of ^ 10^ yr 
(jKrausdliooi ). IPrivon et all (|2012l ) derive an overah SFR 
of 10 to TOMoyr-^ for Cy gnus A. The morphology of the 
young stars is strikingly similar to the one found in our sim- 
ulations. The radius of the stellar ring of ~ 2 kpc would 
mean that the ISM is rather extended in the vertical direc- 
tion, but there are currently no observations available con- 
straining the details of the distribution of a dense gaseous 
component. 

Despite the small number of examples, the action of 
positive feedback is well compatible with or even supported 
by a number of observations mentioned in Sect. [J In partic- 
ular, the recent star formation activity i n more than 75 per 
cent of small (< 15 kpc) radio galaxies (jPicken et al.ll20l"3 l 
or the young stellar populations found in radio galaxies 
(iTadhunter et al.1 12002: IWills et al.ll2002l: [Holt et al.ll2007l : 



iBaldi k CapettiH2008l : jTadhunter et al. ''20H") and the cen- 
tres of radio galaxies ([Aretxaga et al.ir2001 ) may very well 
be caused by jet-induced feedback. And although the young 
stellar populations are mostly discussed in the context of 
merger-induced jet activity, it is not clear whether the 
merger scenario is a viable scenario and could explain the 
observations satisfactorily, whether it can be distinguished 
against jet-induced star formation and what effect jet activ- 
ity would have while a merger is ongoing. The mere fact that 
the merger, increased star formation and radio jets occur at 
similar times does not give a causal relation, in particular 
since the duration of a single jet activity cycle is on the or- 
der of or even below the age accuracy of the young stellar 
population. 

The additional stars in our simulation that have 
been formed due to the action of the jet are almost 
entirely located within the gaseous galactic disc. This 
may however not be interpreted as evidence against the 
possible jet induced star formation awa y from the disc, 
align e d with the axis o f the radio jet (e.g. iMcCarthv et al.l 
Il987l : iDev et al.1 Il997l ). This possible mode of star for- 
mation is associat ed with tens of k pc sized emission 
line haloes (e.g. Meisenheimer Hippelein 1 992: McCarthvl 
Il993l : iBest et aLl7l997l : iMilev De Breuckl |200^), which 
are quite possibly due to gas extracted from the galaxy 
(iNesvadba et al .11 20081 : fKrause Gaiblerll2010l ). In our simu- 
lation, we find that approximately 10^ Mq of disc gas leaves 
the galaxy during the time the jet is active, though this gas 
does not remain dense enough to form the observed emission 
line gas. The absolute number is rather encouraging, and 
should be regarded as a lower limit, since we neglect several 
important effects in this regard in our simulations: First, we 
do not simulate the initial phase of the jet until it reaches the 
extent of 1 kpc. This ph ase might result in substa ntial out- 
ward motions of the gas (jWagner BicknellllioTTI ). Second, 
the ISM is expected to be turbulent, which might also help to 
get it started from the disc. Third, the enhanced supernova 



feedback due to the jet-induced star formation within the 
disc should also add a considerable amount of kinetic energy 
to the gas. Once having left the disc, the gas would be fur- 
ther stirred up by interacti on with the turbulent jet cocoon 
(jKrause Alexanderll2007l V Also magnetic fields (neglected 
in our simulation) might help to stabilize the emission line 
phase. Finally, while generally speaki ng we simulate a high 
power (FR II, iFanaroff Rile v 1974V the jet power m our 
simulation is set rather conservatively. If we would increase 
it by a factor of ten, and if the amount of ejected mass 
would scale linearly with the jet power, we would reach the 
right or der of magnitude that h as been claimed by the obser- 
vations (jNesvadba et al.|[2008l V Since there is not even gas 
cold enough to form emission line regions present in our sim- 
ulations outside the galactic disc, which is however clearly 
present in nature, it is evident, that our simulations do not 
address jet-induced star formation outside the galactic disc. 

Our simulations include a high-power FR II jet - a case 
which is quite rare in the local Universe. The low-power (FR 
I) jets are much more abundant and evolve differently. Yet, 
the early stage may be very similar since FR I jets are still 
relativistic at small scales and would deposit their mechani- 
cal power into the ISM in a similar way as for the simulated 
case. Since their luminosities are much lower, their effect 
on the dense gaseous phase might be much smaller and re- 
stricted rather to small radii, which is possibly difficult to 
detect. In general, low-redshift radio galaxies show low star 
formation rates as derived from poly cyclic arornatic hydro- 
carbo n (PAH) emission features (Ishi et al.ll2007l : IOgle etajj 
I2OIOI ). even though some of them have fairly large masses of 
molecular gas. Yet, we believe that a careful examination of 
radio galaxies for a central star formation and a global en- 
hancement of star formation once the bow shock has prop- 
agated beyond the galaxy should be done, including forma- 
tion of molecular gas at these locations, in order to shed 
more light on the impact of jet activity on the evolution of 
their host galaxie s. Observations at h igh redshift are clearly 
more demanding. IStacev et al ] (|2010l ) detected 158/im [CII] 
line emission in galaxies at 2; 1 — 2, and inferred most of the 
FIR and [CII] emission originates from kpc-scale cores in the 
AGN-dominated galaxies where the local FUV field is much 
higher than in the star-formation-dominated galaxies. This 
led them to assert that the onset of AGN activity stimulates 
large-scale star formation, an argument also consistent with 
detection of strong 88 /xm [OI II] and 122 jtxm [NH] emissio n 
in composite starburs t AGN (|Ferkinhoff et al.ll2010l . l201lh . 
ISevmour et"all (I2OO8I ) were able to detect the PAH signa- 
ture of a large star formation rate in a high-redshift radio 
galaxy. The wide range of observational results found for 
radio galaxies sources might indicate that positive feedback 
depends on certain conditions, and possibly be more effec- 
tive for large disc masses which are more often found at 
higher redshift. 
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